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Abstract 

The reaction 77 J/^J/^ is discussed assuming dominance of the QCD 
BFKL pomeron exchange. We give prediction for the cross-section of this process 
for LEP2 and TESLA energies. We solve the BFKL equation in the non-forward 
configuration taking into account dominant non-leading effects which come from 
the requirement that the virtuality of the exchanged gluons along the gluon ladder 
is controlled by their transverse momentum squared. We compare our results 
with those corresponding to the simple two gluon exchange mechanism and with 
the BFKL pomeron exchange in the leading logarithmic approximation. The 
BFKL effects are found to generate a steeper t-dependence than the two gluon 
exchange. The cross-section is found to increase with increasing CM energy W 
as {W'^)'^^. The parameter A is slowly varying with W and takes the values 
A ~ 0.23 — 0.28. The magnitude of the total cross-section for the process 77 — > 
J/^J/^ is found to increase from 4 to 26 pb within the energy range accessible 
at LEP2. The magnitude of the total cross-section for the process e~^e~ — > 
e^e^J/^'J/^' with antitagged and e~ is estimated to be around 0.1 pb at 
LEP2. 
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The high energy hmit of elementary processes in perturbative QCD is at present 
theoretically fairly well understood [|I|, 0. The leading behaviour is controlled by the 
pomeron singularity which corresponds to the sum of ladder diagrams with reggeized 
gluons along the chain. This sum is described by the Balitzkij, Fadin, Kuraev, Lipa- 
tov (BFKL) equation 0]. Possible phenomenological tests of the perturbative QCD 
pomeron exchange are however difficult. First of all they have to be limited to (semi-) 
hard processes in which the presence of the hard scale (s) can justify the use of pertur- 
bative QCD. Moreover in order to minimize the possible role of the non-perturbative 
contributions it is in principle necessary to focus on the processes which directly probe 
the high energy limit of partonic amplitudes alone. Finally in order to extract the 
genuine BFKL effects which go beyond the conventional QCD evolution with ordered 
transverse momenta from one scale to another it is also useful to consider the processes 
with small (or equal to zero) "evolution length" i.e. those where the magnitudes of the 
two hard scales are comparable. The two classical processes which can probe the QCD 
pomeron by fulfilling these criteria are deep inelastic events accompanied by an ener- 
getic (forward) jet [||, |^) and the production of large pt jets separated by the rapidity 
gap 0. The former process probes the QCD pomeron in the forward direction while 
the latter refiects the elastic scattering of partons via the QCD pomeron exchange with 
non-zero (and large) momentum transfer. Another possible probe of the QCD pomeron 
at (large) momentum transfers can be provided by the diffractive vector meson photo- 
production accompanied by proton dissociation in order to avoid nucleon form-factor 
effects [0, H, while the complementary measurement to deep inelastic scattering + jet 
events may also be the total 7*7* cross section of virtual photons having comparable 
virtuality 0. 



In this paper we wish to analyze the double diffractive production of J/\E' in 77 
collisions i.e. the process 77 — J/\E'J/\E' assuming exchange of the QCD pomeron (see 
Fig. 1). It should be noted that both sides of the diagram shown in Fig. 1 are charac- 
terized by the same (hard) scale provided in this case by the relatively large charmed 
quark mass. In this sense this process is complementary to the classical measurements 
listed above. One of its merits is the fact that in this process we can in principle "scan" 
the perturbative QCD pomeron for arbitrary values of the momentum transfer. In fact 
the diffractive reaction 77 J/\E'J/^ is unique in this respect since the measurements 
listed above do only probe the QCD pomeron either in forward direction or for large 
momentum transfers. In the diffractive double J/^ production in 77 collisions the 
hard scale is provided by the charmed quark mass and the theoretical description in 
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Figure 1: The QCD pomeron exchange mechanism of the process 77 J/^J/^. 

terms of the perturbative QCD pomeron exchange is appUcable for arbitrary momen- 
tum transfers. This process has also the advantage that its cross-section can be almost 
entirely calculated perturbatively. The only non-perturbative element is a parameter 
determined by the J/\Ef light cone wave function which can however be obtained from 
the measurement of the leptonic width rj/5,_»;+;- of the J/\l/. 



The imaginary part ImA{W'^ ,t = —Q"^) of the amplitude for the process 77 — > 
J/'^J/'^ which corresponds to the diagram in Fig. 1 illustrating the QCD pomeron 
exchange can be written in the following form: 

<fk %{k\Q^)<!>{x,k,Q) , , 



ImA{W\t = = j 



n [(fc + Q/2)2 + so][(fc-Q/2)2 + so] 

In this equation x = m'^/i^/W'^ where W denotes the total CM energy of the 77 sys- 
tem, mj/ij, is the mass of the J/\E' meson, Q/2 ± k denote the transverse momenta of 
the exchanged gluons and Q is the transverse part of the momentum transfer. In the 
propagators corresponding to the exchanged gluons we include the parameter sq which 
can be viewed upon as the effective representation of the inverse of the colour confine- 
ment radius squared. Sensitivity of the cross-section to its magnitude can serve as an 
estimate of the sensitivity of the results to the contribution coming from the infrared 
region. It should be noted that formula (|l]) gives finite result in the limit sq = 0. 



2 






Figure 2: The diagrams describing the couphng of two gluons to the 7 — > J/^' transition 
vertex. 

The impact factor $o(^^; Q"^) describes the 7J/\1/ transition induced by two gluons 
and the diagrams defining this factor are illustrated in Fig. 2. In the nonrelativistic 
approximation they give the following formula for $o(^^, Q"^) [05 IHl- 



(J 



1 



1 



g2 m^/^/4 + A;2 



where 



C = qc-^'n-mj/q,fj/^ 



with Qc = 2/3 denoting the charge of a charm quark and 



-2 _ ""^J/* ~^ 



(2) 



(3) 



(4) 



The parameter which characterizes the light cone wave function of the J/\E' can 
be related in the leading order to the leptonic width Tj/^^i+i- of the J/"^ 



fj/-f 



I3mj/if,rj/^^i+i- 



(5) 



In our calculations we will set jjj^ = 0.38 GeV. Equation (|]) can in principle 
acquire higher order corrections which would affect the normalization factor C. It has 
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however been argued in ref. [12] that those corrections should be small provided that 
the light cone wave function is consistently used. The apparently large corrections 
are only present in the non-relativistic potential models where they correspond to the 
effect of "undressing" the constituent quarks ||12| . 



The function k, Q) satisfies the BFKL equation which in the leading ln{l/x) 
approximation has the following form: 



27r2 



k\ 



'2 , „ ' I '2 I „ ^ II a 



k^ + So ^2 + ■50 

k\ 



+ 



where 



and 



fcf + {k' - ky + 2so k'^ + {k' - kf + 2so 

kio 




(6) 



1,2 2 



(7) 



denote the transverse momenta of the gluons. The scale of the QCD coupling as which 
appears in equations (0) and (^) will be set /i^ = k"^ + (5^/4 + m^ where rric denotes the 
mass of the charmed quark. The differential cross-section is related in the following 
way to the amplitude A: 

(8) 



The BFKL equation sums ladder diagrams with (reggeized) gluon exchange 
along the ladder. Its kernel contains therefore the virtual corrections responsible for 
gluon reggeization besides the real gluon emission contribution. The former are given 
by that part of the integral in the right hand side of the equation (P) whose integrand 
is proportional to fc, Q) while the latter corresponds to the remaining part of the 
integral. If in eq. (|l|) one approximates the function by the impact factor 

$o(^^5 Q^) then one gets the two (elementary) gluon exchange contribution to the pro- 
cess 77 J/'^J/'^ |Tl|]. The two gluon exchange mechanism gives the cross-section 
which is independent of energy. A possible increase of the cross-section with energy is 
described by the BFKL effects generated by the solution of equation (|). These effects 
can also significantly affect the t dependence of the cross-section and so the process 
77 J I'^J/'^ might be a useful tool for probing the t-dependence which follows from 
the BFKL equation. Let us recall that in the diffractive photo-production of J/\E' on a 
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proton a possible nontrivial t-dependence generated by the BFKL equation cannot be 



detected due to (non-perturbative) coupling of the two gluon system to a proton |TD 



In order to avoid this effect one may consider the diffractive vector meson photopro- 
duction on a proton accompanied by proton dissociation 0, ||. In this case however 
the momentum transfer has to be large. 

It is known that the BFKL equation can acquire significant non-leading contribu- 
tions |jl3|, |l^, |15| . Although the structure of those corrections is fairly complicated their 
dominant part is rather simple and follows from restricting the integration region in 
the real emission term in equation (I 



For Q = the relevant limitation is ||T6| , [T^ 

(9) 



,x 



X 



It follows from the requirement that the virtuality of the gluons exchanged along the 
chain is dominated by the transverse momentum squared. The constraint @ can be 
shown to exhaust about 70% of the next-to-leading corrections to the QCD pomeron 
intercept |T8|. Generalization of the constraint to the case of a non-forward 
configuration with > takes the following form: 



k'' < {k' + QV4) 



X 



X 



(10) 



Besides the BFKL equation in the leading logarithmic approximation we shall 
therefore also consider the equation which will embody the constraint (|T0]) in order to 
estimate possible effect of the non-leading contributions. 



The corresponding equation which contains constraint ([TD|) in the real emission 
term reads: 



^x,k,Q) = %ik\Q') + 



2 „2\ 3Q;s(/i^) /"i dx' 



dPk' 



27r2 x' J (fc'-fc)2 + so 



,'2 , „ + 7J2 



{k' - kf + So 



Xt + so ' k''^ + SQ {k'^ + so){k'^ + So) _ 
<l>{x', k', Q)Q [{e + QV4)x7x - - 



X 



^2 



fc'/ + (fc' - fc)2 + 2so + {k' - kf + 2so 



Hx',k,Q) 



We solved equations (H) and (|Tl|) numerically setting rric = mj/,j/2, h^qcD = 
0.23 GeV and using the one loop approximation for the QCD coupling as with the 
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number of flavours Nf = 4. Brief summary of the numerical method and of the 
adopted approximations in solving equations (§,|T3) will be given below. Let us recall 
that we used a running coupling with the scale /x^ = fc^ + + m^. The parameter 
So was varied within the range 0.04 GeV^ < So < 0.16 GeV^ It should be noted that 
the solutions of equations |Tl|) and the amplitude (|ID are finite in the limit so = 0. 
This follows from the fact that both impact factors ^Q{k'^,Q^) and ^{x,k,Q) vanish 
for k = ±Q/2 (see equations (||, ||, |ll|))- The results with finite So are however more 
realistic. 



For fixed (i.e. non-running) coupling and for So = equations (|6|,[rTD could in 
principle be solved analytically taking advantage of the conformal invariance of their 
kernels. The numerical method which we adopted is however more flexible and allows 
to analyze equations (|6|.[TT|) in the more realistic case of running and non-zero So. 



In order to solve equations (|6|. [Tl|) we at first expand the function $(a;, k, Q) in the 
(truncated) Fourier series 



N 



$(a;, k,Q) = Y^ ^mix, /c^, Q^)cos(2m0) 





(12) 



where (p denotes the azimuthal angle between the (two-dimensional) vectors k and Q 
and then we discretize the corresponding system of integral equations for the functions 
$m(a;, fc^,Q^) by the Tchebyshev interpolation method. We found that the BFKL 
equation (|^) can be to a very good accuracy approximated by retaining only the term 
corresponding to m = that corresponds to neglecting possible dependence of the func- 
tion $(x, k, Q) upon the azimuthal angle 0. After retaining only the term $o(^) k"^, Q^) 
one gets the following equation for this function: 



3a^(/i 



TT 



2^ "^dx' 



dk'^R{k\ k'\ Q^)%{x', k'\ Q2)_ 
V{k\Q^)^o{x',k\Q^) 



(13) 



where 



R{k\k'\Q^^ 



(QV4 + k'^ + so)2 - Q^k'^ y J{k^ + k'^ + so)2 - Ak^k'^ 
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Figure 3: Energy dependence of the cross-section for the process 77 J/^J/^. The 
two lower curves correspond to the calculations based on equation (|l^) which contains the 
non- leading effects coming from the constraint (0). The continuous line corresponds to 
So = 0.04 GeV'^ and the dashed line to sq = 0.16 GeV^. The two upper curves correspond 
to equation ^ i.e. to the BFKL equation in the leading logarithmic approximation. The 
dashed-dotted line corresponds to sq = 0.04 GeV^ and short dashed line to sq = 0.16 GeV^. 
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+ gV4) 



2U2 



and 



V{k\Q^) = / dk'- 
Jo 





(QV4 + k'^ + soY - Q^k'^ + QV4 + k'^ + sq 

(g2/4 + f,f2 + So) j 



1 



(14) 



{k^ + QV4) 



(gV4 + k'^ + soY - Q^k'^ [ ^(F + k'^ + so)2 - Ak^k'^ 

Q'k^ 

(gV4 + k'^ + So)' - Q^k'^ + QV4 + k'^ + So 



2n rirh pl 



fc2 + Q74 + kQcos{(f)) 



27rJo [k^ + Qy4 + kQcos{(j)) + {l + Xyso] 



(15) 



We assume that similar approximation can be adopted in the solution of eq. (|TT]). 
Equation ([TT|) then takes a similar form to equation ([T3|) with additional constraint 
0[(/c^ + Q'^/A)x'/x — k'"^] imposed on the real emission terms i.e. on those terms in 
eq. (p!3| ) in which the corresponding integrands contain the factors ^o{x', k'"^, Q"^). The 
corresponding equation reads: 



3as(/U^ 



TT 



1 dx' 

X X' 



|)o(x, k\Q'') = ^o{k\Q^) + 
dk'^R{k\ k'\ Q^)%{x', k'\ Q2)0 nj,2 ^ Q^i^yi^ _ ),'2\ _ 



V{k\Q')%{x\e,Q 



(16) 



We based our calculations on the solutions of equations (|T3|, 0). In Fig. 3 we show 
the cross-section for the process 77 —>■ J/'^J/'^ plotted as function of the total CM 
energy W . We show results based on the BFKL equation in the leading logarithmic 
approximation as well as those which include the dominant non-leading effects. The 
calculations were performed for the two values of the parameter sq i.e. So = 0.04 GeV' 
and So = 0.16 GeV'. We have also estimated the total cross-section for the process 
e~^e~ — > e^e~ J/'^J/'^ with antitagged and e~ for the LEP2 energies assuming the 
same cuts as in ref. |19[. We get cTg+e-^e+e- j/*j/*(v^ = VJbGeV) = 0.12 pb and 
ae+e-^e+e-j/^j/^iVs = 175) = 0.09 pb for So = 0.04 GeV^ and so = 0.16 GeV^ re- 
spectively. 



In Fig. 4 we show the t-dependence of the cross-section calculated for sq = 0.10 GeV'. 
We show in this Figure results for two values of the CM energy W (W = 50 GeV and 
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W = 125 GeV) obtained from the solution of the BFKL equation with the non-leading 
effects taken into account (see eq. (|T^)) and confront them with the Born term which 
corresponds to the two (elementary) gluon exchange. The latter is of course indepen- 
dent of the energy W. The values of the energy W were chosen to be in the region 
which may be accessible at LEP2. The following points should be emphasized: 

1. We see from Fig. 3 that the effect of the non- leading contributions is very impor- 
tant and that they significantly reduce magnitude of the cross-section and slow 
down its increase with increasing CM energy W. 

2. The magnitude of the cross-section decreases with increasing magnitude of the 
parameter Sq which controls the contribution coming from the infrared region. 



This effect is however much weaker than that generated by the constraint ([lOD 
which gives the dominant non-leading contribution. The energy dependence of 
the cross-section is practically unaffected by the parameter sq. 

3. It can be seen from Fig. 3 that the cross-section exhibits an approximate (VT^)^^ 
dependence. The parameter A slowly increases with increasing energy W and 
changes from A ^ 0.23 at = 20 GeV to A ^ 0.28 at W = 500 GeV i.e. 
within the energy range which is relevant for LEP2 and for possible TESLA 
measurements. These results correspond to the solution of the BFKL equation 
( p!6D which contains the non-leading effects generated by the constraint (0). The 
(predicted) energy dependence of the cross-section {{W'^)'^^,\ ~ 0.23 — 0.28) is 
marginally steeper than that observed in J/\l' photo-production It should 



however be remembered that the non-leading effects which we have taken into 
account although being the dominant ones still do not exhaust all next-to-leading 
QCD corrections to the BFKL kernel [|13|. The remaining contributions are 
expected to reduce the parameter A but their effect may be expected to be less 



important than that generated by the constraint ([T0|) . The BFKL equation in the 
leading logarithmic approximation generates a much stronger energy dependence 
of the cross-section (see Fig. 3). 

4. The enhancement of the cross-section is still appreciable after including the dom- 
inant non-leading contribution which follows from the constraint (p!0|). Thus 
while in the Born approximation (i.e. for the elementary two gluon exchange 
which gives an energy independent cross-section) we get atot ~ 1-9 — 2.6 pb the 
cross-section calculated from the solution of the BFKL equation with the non- 
leading effects taken into account can reach the value 4 pb at = 20 GeV and 
26 pb for W = 100 GeV i.e. for energies which can be accessible at LEP2. 
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5. Plots shown in Fig. 4 show that the BFKL effects significantly affect the t- 
dependence of the differential cross-section leading to steeper t-dependence than 
that generated by the Born term. Possible energy dependence of the diffractive 
slope is found to be very weak (see Fig. 4). A similar result was also found in 
the BFKL equation in the leading logarithmic approximation [||. 

6. The variation of the parameter sq within the range 0.04 GeV^ < sq < 0.16 GeV^ 
changes normalization of the cross-section by about 30 %. There may still be 
other sources of normalization uncertainties coming for instance from the use of 
the nonrelativistic approximation of the impact factor etc. which can increase the 
normalization error up to 50% or so. The energy dependence of the cross-section 
is however an unambigous theoretical prediction. 

In our calculations we have assumed dominance of the imaginary part of the pro- 
duction amplitude. The effect of the real part can be taken into account by multiplying 
the cross-section by the correction factor 1 + tg'^{Ti\/2) which for A ~ 0.25 can intro- 
duce additional enhancement of about 20 %. 



It may finally be instructive to confront our results with recent findings concerning 



the solution of the BFKL equation in the next-to- leading approximation [0, 14, 15 
It has been found that the non leading effects are very important and that the effec- 
tive intercept A can become negative in the next-to-leading approximation for relevant 
values of as > 0.15. It has been argued that the next-to-leading approximation is not 
reliable and that one has to perform complete resummation of the non-leading contri- 
butions 0, 0. Let us observe that equation (|16D resums to all orders non- leading 
effects generated by the contraint ([ToD . The difference between exact solution of this 
equation and its next-to-leading approximation was discussed in ref. [0 for t = 
and for the fixed coupling a^. It has in particular been found that the exponent A 
corresponding to the exact solution stays always positive for arbitrary values of the 
coupling as- The next-to- leading approximation for this exponent differs significantly 
from the exact solution already for as > 0.2 and can again become negative. This 
result confirms the observation 0, |15| that the next-to-leading approximation alone is 



unreliable and that one has to perform complete resummation of the non-leading effects. 



To sum up we have developed the formalism that enabled us to estimate the con- 
tribution of the QCD pomeron in the 77 J/'^J/'^ diffractive production process. 
We found that the BFKL effects give a significant enhancement of the cross-section 
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Figure 4: The differential cross-section of tlie process 77 J/^J/^ corresponding to the 
solution of equation (|^) which contains the non-leading effects coming from the kinematical 
constraint (|lO[) shown for two values of the CM energy W,W = 50 GeV (continuous line) and 
W = 125 GeV (dashed line). The short dashed line corresponds to the Born term i.e. to the 
elementary two gluon exchange mechanism which gives the energy independent cross-section. 
The parameter sq was set equal to 0.10 GeV^. 



11 



and modify the t-dependence of the Bom term. The cross-section exhibits an ap- 
proximate power law dependence {W"^)"^^ with A ~ 0.25. We based our calculations 
on the BFKL equation which contained non-leading contributions coming from the 
constraint imposed upon the available phase space. This constraint is the dominant 
non-leading effect and it exhausts about 70% of the next-to-leading corrections to the 
BFKL pomeron intercept. We found that the non-leading contributions generated by 
constraint ([lOD significantly affect theoretical expectations based on the BFKL equa- 
tion in the leading logarithmic approximation. This means that the enhancement of 
the cross-section although still quite appreciable should be much smaller than that 
which follows from estimates based on the leading logarithmic approximation |19]. 
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